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Modern weather and climate prediction models are based on a system of nonlinear par¬ 
tial differential equations called the primitive equations. Lie symmetries of the primitive 
equations are computed and the structure of its maximal Lie invariance algebra, which is 
infinite dimensional, is studied. The maximal Lie invariance algebra for the case of a nonzero 
constant Coriolis parameter is mapped to the case of vanishing Coriolis force. The same 
mapping allows one to transform the constantly rotating primitive equations to the equa¬ 
tions in a resting reference frame. This mapping is used to obtain exact solutions for the 
rotating case from exact solutions from the nonrotating equations. Another important re¬ 
sult of the paper is the computation of the complete point symmetry group of the primitive 
equations using the algebraic method. 


1 Introduction 

A main motivation for the study of Lie symmetries of differential equations is that they provide 
systematic tools which allow finding of ansatzes that reduce the number of independent variables 
in partial differential equations. Depending on the particular form of the reduction ansatz, the 
reduced differential equations can then be often integrated to yield exact solutions that are also 
particular solutions of the initial system of partial differential equations. By their definition, 
point symmetries also can be used for generating new exact solutions from known ones. 

Another important application of symmetries of differential equations is that they can pro¬ 
vide a necessary condition of whether two equations can be mapped to each other. Thus, this 
criterion is effective in the case when the target equation is linear as then the initial equation 
is linearizable. For the equations of hydro-thermodynamics, Lie symmetries have proved to be 
extremely successful in finding point, contact and even nonlocal transformations relating dif¬ 
ferent equations. A number of equations, such as the one-dimensional shallow-water equations, 
the Thomas equation and the potential Burgers equation, are linearizable by point transforma¬ 
tions [SlEllIZKlg]. Some equations are linearizable in a nonlocal way, e.g., by a point transfor¬ 
mation after introducing potentials. The most famous example for a nonlocal transformation 
is certainly the linearization of the Burgers equation by means of the Hopf-Cole transforma¬ 
tion [2D], which in fact is a non-invertible transformation, mapping the Burgers equation to the 
linear heat equation. The (—2)-power diffusion equation and some nonlinear wave equations are 
also linearized by point transformations after potentialization. Nonlinear differential equations 
can also be reduced by point transformations to other nonlinear differential equations of sim¬ 
pler form, e.g., the cylindrical Korteweg-de Vries equation to the classical Korteweg-de Vries 
equation. The Liouville equation is linearized by a differential substitution. All the above trans¬ 
formations can be found by invoking the structure of the maximal Lie invariance algebras of 
the equations involved. For invertible point transformations as will be considered in the present 
paper, the relevant necessary criterion for the existence of a mapping relating two system of 
differential equations to each other is that the maximal Lie invariance algebras of the initial and 
the target system are isomorphic 13 Ej. 
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Quite recently, point transformations allowing for canceling terms related to the Coriolis 
force were found for a number of models of fluid dynamics. Although somehow expectable from 
the physical point of view, these transformations are often nontrivial. Examples of particular 
models where such a transformation was already found are the vorticity equation in spherical 
coordinates m n [22], the barotropic potential vorticity equation [2] and the shallow-water 
equations on flat [TO] and parabolic topography Ha¬ 
lt is the purpose of the present paper to show that a transformation eliminating the Coriolis 
force also exists for the more complex system of the primitive equations, which are nonlinear 
partial differential equations for the momentum, mass and energy conservation of atmospheric 
flows. The primitive equations form the dynamical core of most of the modern large-scale 
weather and climate prediction models. 

A further major result of the present paper is the computation of the complete point symme¬ 
try group of the primitive equations using the version of the algebraic method proposed in |3] 
and further developed in m- The method applied can be seen as a refinement of the general 
algebraic technique suggested in imiiT] (see also m) by using the notion of megaideals (i.e., 
characteristically nilpotent ideals) of Lie algebras [25]. To the best of our knowledge, this is the 
hrst example of computing, within the framework of the algebraic approach, the complete point 
symmetry group for a multidimensional system of nonlinear partial differential equations whose 
maximal Lie invariance algebra is inhnite-dimensional and is of complicated structure. The core 
part of the computation procedure is the construction of a sufficiently wide set of megaideals that 
place suitable restrictions on admitted point symmetries of the primitive equations. Without 
the use of megaideals, hnding the complete point symmetry group would require the solution of 
a cumbersome nonlinear system of partial differential equations, which in general is a hopeless 
endeavor. 

The further organization of this paper is as follows. In Section [2] the primitive equations 
are introduced. In Section [3] we compute the Lie symmetries of the primitive equations and 
explicitly find a point transformation that allows canceling of the effects of a constant rotation. 
In Section m we determine the complete point symmetry group of the primitive equations using 
the algebraic method. Section [5| is devoted to the usage of the transformation found and the 
computation of selected exact solutions of the primitive equations. The final Section [6] briefly 
sums up the results of the paper. 

2 The primitive equations 

We consider the primitive equations on the plane using pressure coordinates, i.e. the pressure 
p is used as a vertical coordinate instead of the actual geometric height z. The advantage of 
using pressure coordinates is that the continuity equation reduces to a diagnostic equation in 
this case. The system of primitive equation then reads [18] 

vt + v ■ Vv -I- ujVp + f{-v, u)'^ + V4> = 0, 

‘t'p H- T = 0, 

p 

_n 

“ 1 “ — 0 ; 

Tt + V VT + (uTp - --T = -, 

Cp p Cp 

where v = {u,v) is the horizontal component of the velocity vector, V = {dx,dy) is the two- 
dimensional nabla operator, ui is the vertical velocity in the pressure coordinate system, i.e. 
the material derivative of the pressure p, 4> is the geopotential and T is the temperature. All 
the unknown functions, v, ui, cj) and T depend on {t,x,y,p). Subscriptions of functions denote 
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differentiation with respect to the corresponding variables. The constants /, i?, Cp in the above 
system are the Coriolis parameter, the gas constant for dry air and the specific heat of dry air 
at constant pressure. The function J = J{t,x,y,p) is the external heating. The first equation is 
the momentum equation, the second equation is the hydrostatic equation, the third equation is 
the continuity equation and the last equation is the a version of the first law of hydrodynamics. 
From the physical point of view, the system ([1]) forms the dynamical core of most of the present 
day’s atmospheric numerical models. 

The physical constants R and Cp are always positive. Moreover, from the practical point of 
view, the thermodynamic relation Cp = + R applies for ideal gases. By definition, c„ is the 

specific heat at constant volume, which is the amount of energy needed to heat one kilogram 
of a compound by one Kelvin while holding the volume constant. As there is no compound 
which will be heated by one Kelvin without supplying energy (i.e. always Cn >0), this implies 
that Cp > R. So as to simplify the subsequent expressions, we will put k = R/cp subsequently 
where 0 < k < 1. For the Coriolis parameter / we distinguish between the cases of / = 0 
(no rotation of the reference frame) and / = const (constant rotation of the reference frame). 
There arises the question of whether the choice / = const is a physically interesting one. By 
definition, / = 20 sin y?, where O is the angular velocity of the Earth and p is the geographic 
latitude. Thus, / = f{p) and therefore changes along the meridians. On the other hand, this 
change is rather small and eventually can be neglected for domains extending only moderately 
in North-South direction. For example, for a domain extending approximately 300 kilometer in 
North-South direction, the relative change in the value of / from the South to the North is only 
about 5% around the mid-latitudes. Therefore, for processes that take place on relatively small 
domains, / = const to a good approximation. 

One process that can be described with the model ([T]) for / = const is the land-sea breeze. 
This is a circulation often induced by differential heating of a land-sea boundary, with winds 
directed landward during day and seaward during night. As the land-sea breeze can persist for 
several hours, the effect of the Coriolis force cannot be neglected. It is generally found that 
around six hours after the beginning of the sea breeze the circulation is weakened due to the 
effects of the Coriolis force [ 12 ] . This is why / = const is essential in numerical models that aim 
to capture the land-sea circulation in an accurate way, see [El ED and references therein for a 
detailed review over numerical studies of this particular circulation pattern. 

Another reason why it convenient to assume / = const in the above system is the usage 
of Cartesian coordinates. For processes taking place on a large enough domain, the tangential 
plane approximation of the Earth is not reasonable any more. Eor such processes or for the 
general description of the global atmospheric circulation, it is more appropriate to study the 
primitive equations in spherical coordinates and to use / = 212 sin (p without approximation. 

Within the framework of group analysis of differential equations, the parameterized system © 
should be interpreted as a class of systems of differential equations with the arbitrary elements 
/, R, Cp and J. Two of the arbitrary elements are inessential, in that it is possible to scale 
R = 1 (by a scaling of T and Cp) and / = 1 if / / 0 (by a scaling of {t,x,y,p)). Eor physical 
reasons we will not make a use of these scaling. The possibility to set / = 1 is also not overly 
relevant as we will show in Section [3] that / can be set to zero by a point transformation. 

In what follows we will mostly be concerned with the system © in the case of J = 0 
corresponding to a non-heated atmosphere (the adiabatic case). 

3 Lie symmetries 

We now compute Lie symmetries for the adiabatic case when J = 0 using the infinitesimal 
invariance criterion [20|. We look for infinitesimal generators of one-parameter local symmetry 
groups of the primitive equations ([T|), which constitute a Lie algebra 0 / called the maximal Lie 
invariance algebra of these equations. The subscript / indicates that in fact the algebra Qf 


3 


depends on the the Coriolis parameter /. Each of the infinitesimal generators is a vector field 
of the form 

Q = Tdt + + edp + r]^du + + rj^dr, 

whose coefficients satisfy the system of determining equations implied by the infinitesimal in¬ 
variance criterion. We have computed the algebra 0 / using the Maple package DESOLV [711^177]. 
Splitting the general expression for Q with respect to parameters gives the following basis ele¬ 
ments of 0 /: 

Vi = tdt -f ftydx - ftxdy - {u- ftv - fy)du - {v + ftu + fx)dx - uid^j 
-i2(P + f^{x^ + y^))d^-2TdT, 

V2 = xdx + ydy + udu + vdx + 2(j)dfj) + 2TdT, = pdp + ( 2 ) 

J =-ydx + xdy-vdu + udx, V = dt, S = (cpd^ - Ot), 

^”(7) = 7 • Sx + - htt ■ X + /(7t y - ltx))d4>, Z{a) = ad^, 

where / = // 2 , x = {x,y), = {dx,dy), = {du,dy), 7 = ( 7 ^ 7 ^), and the parameters 

7 ^ = 7 ^(t), 7 ^ = 7 ^(t) and a = a{t) run through the set of smooth functions depending on t. 

The associated one-parametric groups are (i) time-translations, (ii)-(iv) scalings, (v) planar 
rotations, (vi)-(vii) generalized Galilean transformations, (viii) gauging of the geopotential and 
(ix) generalized shifts. 

The algebra 0 / is not singular in the Coriolis parameter /, which means that it is possible to 
set / = 0 in ([2]) . The remaining question is whether there are additional inhnitesimal generators 
extending the algebra 0 o when / = 0 in ([T]). Computing Lie symmetries of system ([T]) for / = 0 
shows that this is not the case, i.e. the primitive equations in a nonrotating reference frame 
admit 0 o as the maximal Lie invariance algebra, whose basis elements are 

Vi = tdt — udu — vdy — uiduj — 24)8^ — 2T St, 

'D2 = xdx + ydy + udu + vdy + 24>d^ + 2TdT, V3 = pdp + 

V = dt, J =-ydx + xdy-vdu + udy, S = p'^{cpd^ - Ot), 

T(7) = 7 • + 7t • - 7it • ^9^, = Oid^, 

The nonzero commutation relations among basis elements of 0 o are exhausted by 

[Dudt] = -dt, [DuS] = 2S, [Pi,T( 7 )] =T(t 7 j, [Vi, Z{a)] = Z{2a + tat), 

[D 2 ,S] = -2S, [P 2 ,T( 7 )] = -T( 7 ), [V2,Z{a)]=-2Z{a), [P3,5] = a^5, 

[dt, T(7)] = T(7t), [dt,Z{a)] = Z{at), [J, X{l)] = - 7 '), 

[T(7), T(cr)] = Z{a "^tt-l' 

Based on the above commutation relations, one can see that the Lie algebra 0 o has the structure 
of 00 = ( 02 © 03 )®i, where 02 = {dt,'Di) is a realization of the two-dimensional nonabelian algebra, 
03 = { 1 ) 2 , V^, J) is a realization of the three-dimensional abelian algebra and \ = {X ( 7 ), Z(a), <S) 
is an infinite dimensional ideal in 0 o, i = ((^’( 7 )) S (Z(q;)))©(5), and {Z{a)) and (5) are abelian 
ideals in the entire algebra 0 o. 

Upon redefining the basis elements in the algebra Qf according to 
dt^dt- fj, T(7) ^ T(7) 

where 7 = ( 7 ^ cos(/f) — 7 ^ sin(/t), 7 ^ sin(/f) -|- 7 ^cos(/t)) and the remaining basis elements 
remain unchanged, they satisfy the same commutation relations as the basis elements of 0 o. 
Therefore, the algebras 0 / and 0 o are isomorphic, which is a necessary condition for the existence 
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of a point transformation mapping the primitive equations with / 7 ^ 0 to the primitive equation 
in a resting reference frame (/ = 0) [5]. This allows us to use the algebraic method for finding 
the transformation relating the two systems with / = 0 and / 7 ^ 0 to each other. 

Suppose that a point transformation 

T: = T{t,x,y,p,u,v,u},(t),T), (4) 


where i € {t,x,y,p,u,v,uj, 

{z\z^,zy,zP,z^,z\z‘^,z^,z'^) 


{t,x,y,p,u,v,uj,(j),T), 

{i,x,y,p,u,v,u,^,f), 


realize the above automorphism between the algebras Qf and 0o- Then, the relation upon the 
corresponding basis elements Q and Q of the algebras Qj and 0 O) respectively, reads 

QV = Qz\ (5) 


which is the usual rule for the transformation of vector helds. 

Evaluating Q for the transformation component of p, it follows from the transformation of Z 
that = 0 and from S that = 0. As T’( 7 )T^ = 0 must hold for arbitrary smooth functions 
7 , one obtains that 7^ = 7^ = 7lf = 7lf = 0. Moreover dtT^ = 0 and from the two different 
scaling operators Vi and P 3 we conclude that p = p- Using similar arguments, one also finds 
that uj = oj, T = T and i =t. 

Now consider the transformation components of x and y. In view of the operators Z and S 
we again conclude that = Trf = = 0. From the condition that = 

and X{'y)T'^ = it follows that T* = a^{t)x + a^{t)y + a^{t)u + a^{t)v + x^{p,uj) and 
= b^{t)x + b‘^{t)y + b^{t)u + b‘^{t)v + y^{p,uj). The condition that {dt — fJ)T^ = 0 implies 
that al — fa^ = 0, aj + fa^ = 0, — /a^ = 0 and + fa? = 0. From JT^ = —T^ it 

follows that b^ = —a^, b^ = a^, b^ = —a'^ and b‘^ = a^. The equation ViT^ = 0 enforces 
that = 0 and therefore also 6 ^ = 6 ^ = 0. Integration of the equations for a? and 

gives that = cicos(/f) + C 2 sin(/t) and = C 2 Cos(/t) — cisin(/t). From the constraint 
a^ 7 ^ + 0 ^ 7 ^ = = cos(/t) 7 ^ — sin(/t) 7 ^ we conclude that ci = 1 and C 2 = 0. The action 

of the scaling operators Di and P 3 finally requires that x^ = y^ = 0. Therefore one obtains 
X = cos{ft)x — sm{ft)y and y = sin(/t)x + cos{ft)y. 

We next consider the transformation components of u and v. Again the condition that 
Z(a)r“ = Z{a)T^ = ST^ = ST^ = 0 leads to = Q. The condition that 

A’( 7 )T“ = aj: and X{‘y)T^ = cr| implies that T“ = d^{t)x + d?{t)y + d^{t)u + d'^{t) + u^{p,uj) and 
{t)x+e‘^{t)y+e^{t)u+e^{t)+v^{p, u}). The relations dj — fd? = 0, d'f+fd^ = 0, d^ — fd^ = 
0 and df + fd? = 0 and = d^, = —d^ and = d^ follow from the transformations 

rules {dt — fJ)T^ = 0 and = —T^, respectively. The scaling transformation T>iT“ = —T“ 
enforces the relations the relation d'^ = d?/f and d^ = —d?/f. Integrating the equations for 
di and d 2 , including the constraint that A’( 7 )T“ = o’} and invoking the transformations of the 
scaling D 3 leads to the transformation u = cos{ft)u — sin{ft)v — f{sm{ft)x + cos{ft)y) and 
V = sin{ft)u + cos{ft)v + f{cos{ft)x — sin{ft)y). 

It thus remains to determine the transformation behavior of cj). From Z{a)T^ = d and 
ST'^ = p'^Cp it follows that = 1 and = 0. The actions {dt — fJ)T^ = 0 and JT'^ = 0 
imply that = 0. From X{‘y)T^ = —^tt ' ^ it follows that Tu = Tv = 0 and Tx = f'^x and 
= Py. The action of the scaling operators T>i and P 3 on T'^ implies that Tp = Tif = 0. 
Therefore, ^ = (f + P{x‘^ + y^)/2. 

It can be checked that all the equations from the condition QT^ = QZ not used to derive 
the above form of the transformation Z = T^{t,x,y,p,u,v,uj, (j),T) reduce to identities. It can 
also be checked by direct substitution that the same transformation also relates the primitive 
equations with / 7 ^ 0 to the equations in which / = 0. This proves the following theorem. 
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Theorem 1. The primitive equations © in a reference frame with constant rotation can be 
transformed to the primitive equations in a reference frame at rest, i.e. f = 0 upon using the 
transformation 


i=t, X = cos{ft)x - sm{ft)y, y = sm(ft)x + cos{ft)y, p = p, 

u = cos{ft)u — sin{ft)v — f{sm{ft)x + cos{ft)y), 

V = sm{ft)u + cos{ft)v + f{cos{ft)x — sin{ft)y), 
f2 

Cb = u;, (^ = (/)++ y^), r = T, 


( 6 ) 


where f = //2. The same transformation maps the maximal Lie invariance algebra qj to qq. 

Remark 1. In cylindrical coordinates, {r,9,p), the transformation ([U]) takes the particular 
simple form 

i=t, f = r, 9 = 9 + ^t, p = p, 

f _ f2 _ 

vL = , hf = + ^r, a) = w, ^ = (f>+ T = T, 

2 8 

where and are the velocity components in radial and in azimuthal direction, respectively. 

So as to derive the transformation Q using the algebraic method it was necessary to assume 
J = 0, i.e. the system was required to be adiabatic. This assumption is crucial as for general 
functions J = J{t,x,y,p) the primitive equations ([I|) only admit the gauging operators Z{a) 
and S. The span of these gauging operators is not enough to derive a sufficient number of 
equations for the transformation components ([ 6 ]). On the other hand, one can check the validity 
of this transformation for the case J 7 ^ 0 by direct computation. As the differential operator 
V • V is invariant under the transformation ([ 6 ]) and this is the only term that is transformed in 
the temperature equation, the same transformation also maps the primitive equations for J 7 ^ 0 
in a rotating reference frame to the corresponding system in the resting reference frame with 
possibly another value of J. 

Corollary 1. Transformation ([6]) maps the non-adiabatic (J ^0) system of primitive equations 
in a reference frame with constant rotation to the non-adiabatic system of primitive equations 
in a reference frame at rest with possibly another value of J. 

Remark 2. Owing to the invariance of the advection operator v-V under the transformation (f6|l, 
it is possible to extend the system of primitive equations ([T]) by equations of the form 


S't + V • VS* + LoSp — Q, 

without introducting new nontrivial transformation components for the prognostic variable S 
and the source term Q, i.e. S = S and Q = Q. Examples for physically relevant equations 
of the above form are, e.g., the moisture equation or any equation for a passively transported 
atmospheric tracer. 

Remark 3. In the case Cp = R, the system of primitive equations where we set 11 = 0 
without loss of generality, admits a wider maximal Lie invariance algebra go than go. Additional 
basis elements go in comparison with go are 

TZ{X) = 2\di + Xixdx + Xpydy - 2X^pdp - {Xpu - X^ix)du - {X^v - - 

(dAj-o) + 2Xiip)dcj — ( 2 X 1 ^ + ■^Xiii{x + y ) \ — 2XiTdf, 
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viip) = i^dp + Tp^dcj + —df, 

p 

where A and ip run through the set of smooth functions of t. It is clear that the operator 7 ^(A) 
is a generalization of the usual shifts in f (A = const) and the scaling operator 2Pi + P2 = t)- 
For arbitrary A, TZ{X) can be interpreted as re-parameterization of time. The operator in 

turn is a generalized Galilean boost in p-direction. At the same time, this example is unphysical 
as Cp > R. 

4 Complete point symmetry group 

The above consideration shows that without loss of generality it suffices to carry out group 
analysis of the primitive equations © only for the case / = 0 In this section we find the 
complete point symmetry group of Eqs. ([T]) with / = 0 by the algebraic method proposed 
in [ 3 ] . This method may be treated as an enhancement of the approach suggested in [MllIT] (see 
also | 15 ] ) by embedding the notion of megaideals [ 25 ] , which is a brief name for fully characteristic 
ideals. See |H] for recent advances of the enhanced method. Its main benefit is that it can be 
applied even to systems of differential equations possessing infinite-dimensional Lie invariance 
algebras, which is the case for Eqs. although it can also simplify related computations for 
higher-dimensional Lie invariance algebras. 

The algebra 0 = 0o has the following obvious megaideals: 

0' = (L>, 5 , A( 7 ),^(a)), 0" = (A(7),Z(a)), 0'" = Zg. = {Z{a)), 

Z,, = {S,Z{ 1 )), Zg,n0'" = (^(l)), 

mi = Cg(0") = (I?3,*S,Z(a)), m; = ( 5 ), Cg(mi) = (J, ^(7), ^(a)), 

where a', Zq and C(j(b) denote the derivative and the center of a Lie algebra a and the centralizer 
of a subalgebra b in a, respectively. Here and in what follows the parameters 7^, 7^ and a run 
through the set of smooth functions depending on t. 

To find more megaideals of 0, we apply Proposition I from | 13 ] for various special choices 
of the megaideals io, ii and i2 of 0. This proposition states that the set s of elements from Iq 
whose commutators with arbitrary elements from ii belong to i2 is also a megaideal of 0. Thus, 
for io = 0'", ii = 0' and i2 = Zg/ n 0'" = {Z{ 1 )), we obtain s = {Z{ 1 ), Z{t)) and hence this is a 
megaideal. We reassign the last s as i2 and iterate the procedure with the same io and ii, which 
gives the series of megaideals {Z{ 1 ), Z{t ),... , Z{t^)), n G Nq. 

A convenient choice for io and ii is io = ii = 0 when i2 is varying. Eor i2 = m'^ and i2 = 0" we 
respectively have the megaideals s = {'D^,S) =: m2 and s = {kV2 + 2 V^, J, Z{a)). Then 

Cg'(nx2) = and Cg(m2) = (T>i-bT>2,T>2 + 2T>3, ^”(7), Z(a)) =: m3 are also 

megaideals, as well as Cm3(Zg/ n 0"') = {Vi + 1)2, JjR, ^{l), Z{a)). 

Applying again Proposition 1 from m on the next step, we take io = ii = Cg'(m2) and 
i2 = q'" and derive the megaideal s = (A’(l, 0 ), ^( 0 , l),Z{a)) =: m4. We reassign the last s as 
i2 and iterate the procedure with the same io and ii, which gives the series of megaideals 

(A(I, 0 ), A( 0 , 1 ), A(f, 0 ), A( 0 , f),... , X{e, 0 ), X{ 0 , G), Z{a)), n G No. 

Considering io = 0 and ii = m4 © m'^ with i2 = 0'", we get s = {'Di,V,S,X{'y),Z{a)). 

Some of the above megaideals of 0o can be neglected in the course of computing the complete 
point symmetry group Go of the primitive equations ([T|) with / = 0 by the algebraic method. 
Indeed, the condition G*i C i for a megaideal i may only result in constraints for components 
of point symmetry transformations that are consequences of those obtained in the course of the 
computation with other megaideals. In particular, this is the case if a megaideal i is a sum 
of other megaideals. To optimize the computation, we select a minimal set of megaideals that 
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allow us to easily derive a set of constraints for components of point symmetry transformations 
that is maximal within the algebraic framework. We choose the following megaideals from those 
we have computed: 

(Z(l)), { 2 ( 1 ), Zit)), ( 5 ), {X{l, 0 ),X{ 0 ,l),Zia)), 

0 ), d:’( 0 , 0 ), d:’( 0 , t), A’(l, 0 ), d:’( 0 , 1 ), Z{a)), ( 7 ) 

{J,;Yij),Z{a)), {V,Xh),Z{a)), {V, + J ,V, X{l), Z{a)), 

{V 3 ,S), {Vi,V,S,X{-f),Zia)). 

We additionally ordered the megaideal list in such a way that megaideals heading the list 
give more elementary equations of the form = 0 with some i,j G {t,x,y,p,u,v,u}, (p,T} or 
allows us to specify the expressions for some . 

The general form of point transformations that acts in the space of the independent and 
dependent variables of the primitive equations m is given by Eq. m, where the corresponding 
Jacobian J does not vanish. For a point transformation T to be qualihed as a point symmetry 
of the primitive equations © with / = 0, its counterpart T* push-forwarding vector fields 
should preserve each of the selected megaideals ([ 7 |) of the algebra go. As a result, we obtain the 
conditions 


r*^(l) = = aiZ(l), (8a) 

%Z{t) = = a2Z{i) + 03^(1), (8b) 

%S = p^icpT} - n)d,. = aj, (8c) 

r*T(l, 0) = = X{h% 6°?) + Z(d°i), (8d) 

r*T(0,1) = Vyd.i = T(6?o, 6“) + Z(d02), (8e) 

%X{t, 0) = {tTi + X)d,. = X{b^^i + + bf^) + Z{a^^), (8f) 

r*A(0, t) = {tr^ + T:)d,. = X{b]},i + + 6i0) + (8g) 

%X{t\Q) = ifT^: + 2 tX - 2 xT;)d,. 

= XX^P + bf,i + 620, 622p ^ ^ ^20) ^ -21)^ 

r*T(0,t2) = XX + - 2 y 7 ^) 5 ,-. 

= XXiP + blli + 620, blX + blli + 620) + Z{aX: (8i) 

%J = {xX -yX + - vT:)d,. = a^J + X{X) + Xd^), (8j) 

XP = Xd-z^ = aeV + X{X) + (8k) 

%iPi + V 2 ) = {tX +xX + yX + 0 jX)d,i = + V 2 ) + asP + X{X) + Z{a^), ( 81 ) 

%Vz = {pX + ^X)dzi = «9^’3 + aio^, (8m) 

%V, = {tX - uX - vX - ^X - ‘^<PX - 

= aiiDi J- CL12P + 0135 J- A’('yO) J- ZX^), ( 8 n) 


where i € {t,x,y,p,u,v,uj, ^,T}; arid we assume summation with respect to repeated indices; 
as, s = 6^?, bl^, blj, bf^, 6^| and bl], k,l = 1,2, are constants; 7™ = (7™\7”*^), 

m = 3 ,..., 6, and the parameters dO*, 7™*, 7™* and d™ are smooth functions depending 

on t. 

We will derive constraints on T by sequentially equating the coefficients of vector fields in 
the conditions (l8|) and by taking into account the constraints obtained in previous steps. 

Thus, the condition (l 5 aD directly implies that = oi and 7 ^ = 0 if i 7^ (/>. Then the 
constant ai is nonzero since the Jacobian J does not vanish. The equation ait = 02! + 03 derived 


from the condition (j8bp gives that ^2 7^ 0 and hence the component depends only on t and 
the dependence is affine, 


7^ — Q,ICL2 — Q3®2 


-1 


This completely specifies expression for and also implies that ( 9 ^ = ^a2dt. 

The condition pci) is split into the equations p'^{cpT^ — T^) = P^{cpT^ — T^) = 

a4Cp(T^)^ and CpT^—T'^ = 0 for i ^ 4 >,T. Therefore, = ai{T'^/pY, 7 ^ = Cpai — Cpan{T'P/pY, 
7 ^ = 0 for i 7^ (^, r, and 04 Y 0- 

Considering simultaneously the pairs of the conditions (I8d[) and (l8ep . ( 1 ^ and p^ , as well 
and (IHIl), we derive that 


as 


'~rx _ 
' X 

-bT,, 

v = 

bfv 

7 t = 

dOi, 












7 p = 


Tx = 

V = T 

i T X 






'Y'x _ 
~ 

-b% 

7 ^ = 

^22) 

d02, 


y ’ 





T^ = 

: b\\i 

+ b\l 

-bTit 

, 7 ^ = 

: + ^21 “ 


77 = 

b\\, 

TV 

= bfi, 

n'<t> — 

' U 

d^^ 

T^ = 

- b\li 

+ b\^ 


, V = 

: 622^ + ^22 “ 

^22^5 

f-j-U _ 

'v 

b\l 

JV 
' V 

= ^221 

n-<t> — 

' V 

di2 

t: = 

V = 

-- 0, i - 

= t,p,~ 

w,T; 










bTit^ 

+ “^^iblji -+ 

-b^H- 

bTit) - 

- /)22;2 

+ Okl 

f 

^bij 

t = 2 bll 


k,l = 

: 1,2 


2tQ:^^ — — 2aix = —2b‘^x — 2b‘^y + 

2tQ:^^ — YdP‘^ — 2 aiy = — 2 b^x — 2 b^y + 


The last two equations imply that |6|^| Y 0 (otherwise, the Jacobian J equals zero) and thus 
the transformation components x = T* and y = depend only on {t,x,y). More precisely, in 
terms of the constants b^ we have the representation 

r- = 6??x + ry = 

where ( 3 ^ are smooth functions of t. AsT^ = 7 u = = 0 , we obtain bl]i+bYi —b^^t = 0. 

Then = aia^^b^] and = afa^'^b'^^, i.e., = aia^B^^ and where 

we use the matrix notation B^^ = (6°°), B^^ = (blj) and B^^ = On the other hand, 

— 2 (B^^YB^'^ = — 2 aiE, where E is the 2 x 2 unit matrix, i.e., {B^^^B^^ = afa^^E, which 
implies, e.g., for the ( 1 , l)-entry that (61?)^ + (^12)^ = afa^'^. Therefore, oi > 0 and thus we 
can represent the matrix B^^ in the form 


where O is a 2 x 2 orthogonal matrix. This completes specifying the expressions for T* and T^. 

The representation for B^^ implies 6^5 = 622 ^12 ~ “^21- Using this, we derive from the 

condition ([^ that B^^x = a^x + (7^^, —7^^)''', which gives 05 = 1, / 3 ^(t) = 7^^(t) and / 3 ^(t) = 
—7^^(t). In view of the above equations for derivatives T^i with i,j = u,v, briefly representable 
as = B^^, we also get from the condition ([^ that ai^a2B^^v = a^v + (yp, —Tp)"''- 

Arranging the last equation results to finally specifying the expressions for T“ and , 

= —(6??u + + / 3 pt)), T" = —{b^^u + b^^v + / 3 pt)), 

ai ai 

As the derivatives Tu and Tv may depend only on t, the condition ([^ gives that 7 p = 
Tv = 0 . The variables x and y are involved in the expression of T‘^ only within the summand 

-al 
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The condition (j8kp obviously implies the elementary equations 7 /’ = = 0 and the 

constraint that the transformation component may involve the variable t only via the above 
summand and one more summand that depends only on t. 

Two more elementary equations, TS = = 0 , follows from the condition ([HU). The equation 

implied by ([8l|) for is tT^ = a’j'T^ + ^8) which gives aj = 1 . Then the equation implied for 
takes he form The main feature of T‘^ obtained from ([8l)l is that the term ujTu 

depends only on t, x and y. 

Consider equations yielded by the condition (jHinj). Thus, the equation for is = 
—oiop'^, and hence 7 ^ = 0 . As = a4^{T^/pY, we get that {T^/p)p = 0 , i.e. pTp = T^. The 
equation for is pTp = agT^. Therefore, ag = 1 . Then the equation for 1 '^ takes the form 
pTY + and, after combining with the analogous equation that is obtained from the 


condition (|8l|), reduces to the equation = 0. 

Collecting coefficients of 9 ^, df and in (f8nl) results in the equations tTt = au'Y + au, 
— 2 TTY = — 2 aii 7 ^ — and 

tTt^ - - 4>T^ - 2TT^ = -2aiiT^ + a^CpiTPY - 7?*' x + 

The essential consequence of the first equation is an = 1 . Then the third equation implies that 
does not depend on and oj since we have already proved that all summands in the left hand 
side of the equation as well as have this property. From the second and third equations it 
is obvious he variable p is involved in the expressions of and only within the summands 
and \ai‘iCp{TY'^ 1 respectively. 

We introduce notation of the following constants: 


as / n n A 1 K 

eo =-, ei = — / 0, 62 = >0, 63 = — >0, 64 = - ai 2 el , 

a2 a2 a2 p 2 


ai 


T" 


65 — 7^ 0 , 6g — a 4 S 3 . 

UJ 

The constant 63 should be greater than zero for both physical and mathematical reasons since 
the exponent 63 should be well defined for all k: 0 < «: < 1 and, in view of the physical 
interpretation of the variable p, both its initial and transformed values should simultaneously 
be positive. The constant 62 can be assumed positive since the parameters 62 and O are defined 
up to simultaneously alternating their signs. Collecting all the restrictions we have derived for 
the components of the transformation T within the algebraic approach and using the above 
notation, we obtain the preliminary representation of this transformation. 


i=eit + eo, X = 62OX +/ 3 (t), p = e^p^ 

Cj = 650;, 


V = —Ov + —/3t(t), 

61 61 


( 9 ) 


' + Cjj 


^ - 66 ) T + e^Cpp'^ — ■ Ox + a{t), f = e^T - Eap'^. 


Not all parameters in the representation ([ 9 |) are independent. For the transformation T to 
really be a point symmetry of the primitive equations some parameters have to satisfy 
additional constraints that cannot be derived within the framework of the algebraic approach. 
This is why the computation should be completed by the direct method. The application of the 
direct method can be simplified by factoring out a priori known continuous transformations. 
Thus, we can set 60 = 64 = 0 , /3 = 0, a = 0 and O to be equal to the diagonal matrix with the 
diagonal entries equal to —1 or 1. 

We calculate expressions for transformed derivatives and substitute them to the primitive 
equations © written in terms of the transformed variables, which are with tildes. Then we 
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choose ut, vt, (fip, ojp and Tt as principal derivatives, express them in terms of other (parametric) 
derivatives from ([T|) , substitute the obtained expressions into the system derived on the previous 
step. Splitting the resulting system with respect to parametric derivatives gives the missing 
equations, 

c c-2 

£5 — —! ^6 — -o- 

ei 

This equations jointly with the representation Q leads to the following assertion; 

Theorem 2. The complete point symmetry group of the primitive equations (HD consists of the 
transformations 


i=eit + £o, X = e 20 x +/3(t), p = e^p, 

^ = 1 °'' + “ = 1 “’ ( 10 ) 

^ = -^(p + e 4 Cpp'^ - ■ Ox + a{t), f = -|r - Eip'^, 

El El 

where Eqj • • •, £4 are arbitrary eonstants with ei 7^ 0 , £2 > 0 and £3 > 0 ; (3 = (/ 3 ^,/ 3 ^); the 
parameters and a run through the set of smooth functions of t; O is an arbitrary 2x2 

orthogonal matrix. 

Corollary 2 . The discrete symmetries of the primitive equations (HD are exhausted, up to com¬ 
bining with continuous symmetries and with each other, by two involutions, which are the simul¬ 
taneous inversion of time and veloeity, {t,u,v,w) —)• (—t, —u, —v, —w), and simultaneous mirror 
mappings in the {x,y)- and {u,v)-planes, {x,y,u,v) —>■ {—x,y,—u,v). 


5 Exact solutions 


Finding the transformation ([S]) has two more immediate benefits. It allows one to take arbitrary 
exact solutions of the primitive equations in the resting reference frame to exact solutions of the 
primitive equations in a constantly rotating reference frame and vice versa. This transformation 
is also important because it enables one to carry out Lie reductions using the simplified Lie 
invariance algebra 0 O) spanned by the operators ([ 3 ]) and then to extend the solutions obtained 
to the rotating case. Examples for both of the above usages are presented in this section. 

Physically, the simple solution of the nonrotating primitive equations, 

u = uo{p), V = vo{p), u} = 0, (t) = (pip), T = T{p), 


where the relation between T and p is given via the hydrostatic equation, describes a stably 
stratified atmosphere with a horizontally homogeneous horizontal wind field, a vanishing vertical 
velocity and horizontally homogeneous fields of geopotential and temperature. The inverse of 
transformation ([6]) takes this solution to 


/. 


tt = cos -t uo{p) + sin -t Voip) + -y, 


f. 


f 


f. 


/. 


/ 


u = - sin i^-tj uq{p) + cos \^-t j vq{p) - - x , 
cD = 0, p, = cP{p)-^-^{x^ + y^), f = T{p), 

which is a solution of the primitive equations in a constantly rotating reference frame. This 
solution now is horizontally isotropic in the geopotential, while there is still no vertical velocity. 
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Physically, this means that the effects of a constant rotation cannot lead to vertical motion if 
the initial vertical velocity is vanishing. The above solution then describes the inertia motion 
of fluid particles under the action of the Coriolis force, cf. [TO] for the corresponding solution of 
the rotating shallow-water equations. This type of motion can be frequently observed for buoys 
in the ocean. 

To systematically carry out Lie reductions of the primitive equations ([T|) with / = 0 , it is 
necessary to compute an optimal list of inequivalent subalgebras, which forms the cornerstone of 
the reduction procedure. We do not aim to establish a complete list of inequivalent subalgebras 
of dimensions one, two and three here, which for the proper cases would allow reduction of 
the number of independent variables by one, two or three. In other words, the corresponding 
reduced systems would be systems of partial differential equations in two independent variables, 
systems of ordinary differential equations and systems of algebraic equations, respectively. 

Instead, we consider the Lie reduction with respect to the subalgebra 


s = (T’( 7 ) + aiS,X{a) + a 2 <S), 

where oi and 02 are arbitrary constants, the pairs 7 = (7^,7^) and cr = of smooth 

functions of t are linearly independent and 'y^^-cr — crtf^ = 0 . Note that in this case the operators 
^(7) and X{eT) commute and form a proper subalgebra that is suitable for Lie reduction. 
Previous experience shows, that this algebra is indeed an element of the optimal list of two- 
dimensional subalgebras of the primitive equations, see the corresponding results for the vorticity 
equation, the Euler equations, the Navier-Stokes equations and the magneto-hydrodynamic 
equations [IllllsllMlISS]. An appropriate reduction ansatz corresponding to this subalgebra is 


o"*- • X 7“’- • X 

V = V + —7, - — 

a; = LJ, 

K I I 

^ = $ + Cpy (aicr-L - 027 “^) ' x-' x- 

T = p'^f + y (aicr-^ - 027-^) • X, 

where 7-*- = (7^,-7^), cr^ = a^), 6 = 7^(7^ — 7^17^ = 7 • = —7-*- ■ a 0 {6 can 

be assumed to be positive up to simultaneously alternating signs, e.g., of 7^ and 7^), and 
quantities with hat depend on the invariant independent variables t and p. By the way, for each 
pair P = (/ 3 ^,/ 3 ^) one has the representation 


f 3 


7^-/3 

-r- T- J ^ 




The above ansatz reduces the primitive equations ([T]) with / = 0 to 

~ ~ O’-*- • V 7-*- • V p'^ I I 

vt-bwVpH - - —7i- - —-b Cp-^(aicr -027 ) = 0, 

<^p -b Rp'^~^f = 0, 

— -b CJp — 0, 

ft + ufp + - 027 “^) • V = 0 , 

Upon integrating the reduced continuity equation (|llcp to yield 


(lla) 

(llb) 

(llc) 
(lid) 


w = 6td ^p + x{t), 


12 














it is clear that the above system is reduced to a linear system of four (l+l)-dimensional first-order 
partial differential equations, which can be solved in the following way. We make the change of 
dependent variables v = Gv, where the 2x2 nondegenerate matrix-function G = G{t) is chosen 
as a solution of the equation Gt — HG = 0 with the 2 x 2 matrix-function H = H{t) defined by 
ifx = ( 5 “^(cr-*- • x)')'j — Let G~^ denote the inverse of the matrix G. Then the 

first two equations (lllal) reduce to 

W + wvp-k - 027“^) = 0, (12) 

which is an inhomogeneous system of two decoupled linear partial differential equations. The 
change of the independent variables 

maps the system (I 12 p to the system of trivial ordinary differential equations with the independent 
variable r, where ^ plays the role of parameter. The general solution of the latter system can 
be found by quadratures. This gives the following expression for v: 

v{t,p) = G(f)v °(0 - CpG{t) f (( 5 (r))''"^(C 6l(T))'"G“^(r)(oicr-L(r) - 027-^(r)) dr, 

■Jto 

where is a pair of arbitrary smooth functions of 

In order to solve the equation (jlldp . we substitute the obtained expression for v into it, 
switch again to the variables (r, and integrate with respect to r. As a result, we have 

T{t,p) = r°(0 - ^ ■ v(r,(f(r)(C-h6»(r))) dr, 

where T® is an arbitrary smooth function of 

The last step of solving the system dill) is the integration of the equation dllbp with respect 
to p, which gives 

^ = <(>°(^) -R [ P'^~^T{t,p) dp, 

Jpo 

where is an arbitrary smooth function of t, which can be neglected. 

Substituting the expressions derived for the values with hats into the ansatz, we obtain the 
entire family of s-invariant solutions of the primitive equations ([T]). 

6 Conclusion 

The present paper is devoted to an investigation of the system of primitive equations from its 
symmetry point of view. We found a point transformation that allows to cancel the effects of a 
constant Coriolis force in this system. This transformation might be relevant for application of 
the primitive equations on the /-plane, such as in studies of land-sea breezes. 

In practice, the primitive equations ([I]) as presented in Section [2] are forced and damped 
by external mechanisms, such as external heating (including J), phase transitions of water 
(including an equation for moisture) and bottom friction (including friction in the momentum 
equations). The presence of these additional mechanisms substantially narrows the number of 
admitted Lie symmetries. It was discussed in Section [ 2 ] that for arbitrary J only Z{a) and 
S are admitted as Lie symmetries. Symmetry breaking due to external forcing terms thus 
substantially hinders the applicability of symmetry methods. On the other hand by omitting 
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these external influences it is possible to arrive at a system that has a wide maximal Lie invariance 
(pseudo)group and is therefore accessible to the machinery of group analysis. Eventually, results 
derived for the simplified equations can be extended to the usual system. In this manner we 
have shown in the present paper that the same transformation Q that maps the rotating 
primitive equations to the non-rotating ones in the adiabatic case extends trivially also to the 
non-adiabatic case, without the necessity to modify J in the transformed equation. In other 
words (l6|) in not an equivalence transformation. 

In a similar manner, it can be checked that the transformation Q also maps the dissipative 
primitive equations in a constantly rotating reference frame to the dissipative primitive equations 
in a resting reference frame. This means that attaching classical friction, uAv, to the right-hand 
side of the momentum equations in ([T|) does not require to modify transformation Q in order 
to set / to zero. 

In Section [ 5 ] we discussed an application of the model of primitive equations for the choice 
/ = const, i.e. on the /-plane. For domains extending farther in North-South direction, the 
latitudinal variation of the Coriolis parameter becomes relevant. The next order of accuracy 
approximation for / in a Cartesian plane is / = /o + / 3 y, (5 = const, i.e. a linear variation of 
the rotation. It can be checked that in this case, the primitive equations ([T]) only admit a six- 
parametric maximal Lie invariance group which therefore cannot be isomorphic to the maximal 
Lie invariance algebra go dSl) computed in Section [3l This at once implies that there cannot 
exist a point transformation which maps the case of / = /o -|- / 3 y to the primitive equations in 
a resting reference frame. 

In Section [J] we established another main result of this paper by computing the complete 
point symmetry group of the primitive equations using the algebraic method. This computation 
was rather elaborate due to the wide spaces of both independent and dependent variables of 
the primitive equations. It was necessary to establish a suitable set of megaideals, finding 
of which crucially relied on the iterative use of Proposition 1 from [8]. Without this set of 
megaideals it would have been overly difficult to simplify the determining equations of point 
symmetry transformations for the primitive equations enough to enable their direct integration. 
This example thus shows the power of the algebraic method for finding complete point symmetry 
groups for large systems of nonlinear partial differential equations admitting infinite dimensional 
Lie pseudogroups, which would be challenging with the conventional direct method. 

In Section Owe shortly discussed the construction of exact solutions of the primitive equations 
with rotation. The mapping ([U]) allows one to carry over solutions of the nonrotating equations 
to the equations in a rotating reference frame. This is important from the point of view of exact 
solutions that can be obtained from Lie reduction, as the operators from the algebra go are 
considerable simpler than those from the algebra g/. In the paper, we derived an exact solution 
of the primitive equations that arises from a completely integrable case of Lie reduction. The case 
considered is certainly the most important example of reduction with respect to two-dimensional 
subalgebras. We should also like to stress that upon constructing the optimal lists of inequivalent 
subalgebras, a considerable fraction of these lists will not be suitable for Lie reduction. In 
particular, all algebras including {Z[a)), (S) or some combination of these two basis elements as 
subalgebra will not allow one to And a reduction ansatz. Moreover, a number of cases arising from 
three-dimensional subalgebras will also not be needed for reduction. More precisely, reductions 
using algebras including (^(7) -|- a^S, i^{(T) + a^ 5 ), a^, € M, = 0 for i = 1, 2 as a 

subalgebra are not required the system of differential equations resulting from reduction using 
this subalgebra can be completely integrated by quadratures. 

In view of the remarks of the previous paragraph, despite we have not systematically followed 
the steps of group-invariant reduction, the results obtained are in a certain sense a substantial 
part of the exact solutions of the primitive equations that can be found by Lie reduction. A more 
detailed exposure of the group analysis of the system of primitive equations will be presented 
elsewhere. 
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